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Preface 



These lecture notes have been prepared as a basic introduction to the theory of constrained 
systems which is how basic forces of nature appear in their Hamiltonian formulation. Only 
a basic familiarity of Lagrangian and Hamiltonian formulation of mechanics is assumed. 

The First chapter makes some introductory remarks indicating the context in which various 
types of constrained systems arise. It distinguishes constrained systems for which the equa- 
tions of motion are uniquely specified from those for which the equations of motion are not 
uniquely determined. The focus of the lectures is on the latter types of constraints. The 
notations that will be used is introduced here. 

In the second chapter, the general features are introduced in the familiar example of source- 
free electrodynamics. 

In the third chapter, features seen in electrodynamics are abstracted in the simpler context 
of systems with finitely many degrees of freedom. How constraints arise in a Lagrangian and 
in a Hamiltonian formulation is discussed. 

In the fourth chapter, the Dirac's algorithm for discovering constraints and their classification 
is discussed in a sufficiently general context. 

In the fifth chapter, a special class of constrained systems for which the Hamiltonian itself 
is a constraint is discussed. Such systems arise in the context of general relativity and lead 
to a host of issues of interpretation of 'dynamics'. 

A set of exercises is also included for practice in the last chapter. 

These set of eight lectures were given at the Refresher Course for College Teachers held at 
IMSc during May- June, 2005. 

June 4, 2005 Ghanashyam Date 
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Chapter 1 

Introductory Remarks 



Let us begin by recalling some elementary understanding of what one means by mechanics. 
There are two distinct parts: (a) kinematics and (b) dynamics. Kinematics specifies the 
(possible) states of a given system while dynamics specifies how a system evolves from one 
state to another state. The central problem one wants to solve is to find the state of system 
at some time t given its state at an (earlier) time t . While there are various types of 'state 
spaces' possible, we will be dealing with those systems for which the state space is a manifold 
(finite or infinite dimensional) and an evolution which is of the continuous variety. 

There are three distinct tasks before us (i) devise a framework or a calculational prescription 
which will take as input a state space, some quantity reflecting/encoding a dynamics and 
give us a corresponding rule to evolve any given state of the system. A familiar example 
is the Lagrangian framework wherein the state space is described in terms of positions and 
velocities, dynamics is encoded in terms of a Lagrangian function and the variational principle 
leads to the Euler-Lagrange equations of motion, (ii) Use physical, qualitative analysis of a 
given physical system to associate a particular state space and a particular Lagrangian (say) 
with the physical system and (iii) devise methods to obtain generic evolutions in sufficiently 
explicit terms. 

The last task is where one will make 'predictions' and is the most relevant in applications (eg 
engineering). Here the issues such as sensitivity to initial conditions etc play an important 
role. This task can be addressed only after the first two are specified. This is not the aspect 
we will discuss in these lectures. Some of it will be discussed in the NLDL part of the course. 

To be definite, let us consider the Lagrangian framework. Thus we will have some config- 
uration space, Q, which is some n dimensional manifold with coordinates q\i = 1,2, • • -n. 
The dynamics is encoded in a function L(q\ q\ t) and the equations of motion are the Euler- 
Lagrange equations of motion, 

is = iJ*n q M) = o ^§-§ = 0, v, (i.i) 

Questions: How general can Q be? How general can L be? Can we always obtain equations 
of motion via the variational principle? 

Generally, one takes Q to be an n-dimensional manifold. However, this need not be R N , 
although usually it is. For instance one may begin with n particles moving in 3 dimensions, 
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so that Q is 3n dimensional, but there could be restrictions such as distance between any two 
particles is always fixed (rigid body) or that particles must have velocities tangential to the 
surface of a ball. In the latter case, the relevant Q would be the 2-sphere which is a compact 
manifold without boundary. In both cases the dimension of the relevant configuration space 
is smaller that what we began with. We could also imagine confining the particles to a 3 
dimensional box in which case the relevant configuration space is a bounded portion of R 3n . 
In short, the possible motion of the system may be constrained. A majority of applied me- 
chanics problems have to deal with such constrained systems and there is a massive recent 
treatise on such systems Typically, such constrained systems are described by specifying 
relations f a (t,f,v) = 0. These relations are required to be functionally independent, mu- 
tually consistent, and valid for all possible forces. Some terminology that you might have 
come across classifies the constraints as: (1) Relations independent of velocities are termed 
positional/holonomic/ finite constraints; (2) Velocity dependent ones are called velocity con- 
straints which are further divided into scleronomic (time independent), rheonomic (time 
dependent), holonomic (integrable), non-holonomic (non-integrable). Constrained systems 
of this variety, serve to restrict the configuration space and also possible motions, but the 
Euler-Lagrange equations are always solvable for accelerations, so that dynamics is uniquely 
specified. We will not be dealing with this type of constrained systems. 

The constrained systems that we will deal with arise in situations where the well moti- 
vated choices of the configuration space and the Lagrangian, do not specify the equations 
of motion uniquely. In the Lagrangian framework, these correspond to the so called singu- 
lar Lagrangians. Analysis of such systems is better carried in the Hamiltonian framework 
wherein one arrives at an understanding of gauge theories and we will be essentially focusing 
on such constrained Hamiltonian systems. 

From the point of view of usual applications of classical mechanics, such systems would 
appear quite exotic and possibly 'irrelevant'. However all the four fundamental interactions 
that we know of, when cast in a Lagrangian or Hamiltonian framework, precisely correspond 
to the kind of constrained systems we will discuss. An understanding of these basic forces at 
the classical structural level is crucial also for constructing/understanding the corresponding 
quantum theories. In the context of general relativity, the constrained nature of the theory 
throws up challenging conceptual and interpretational issues. With these reasons as primary 
motivations, we will discuss constrained systems. 

We will also use certain notations which are introduced below. 

1. I will freely use terms such as 'manifold', 'tensors' etc. Here is a very rapid, heuristic 
introduction to these terms. 

The idea of a manifold (actually a different iable manifold), is invented to be able to 
do differential and integral calculus on arbitrarily complicated spaces such as (say) the 
surface of an arbitrarily shaped balloon. The basic definition of differentiation involves 
taking differences of values of a function at two nearby points eg., f(x + h) — f(x), 
dividing it by h and taking the limit h — > 0. This is fine when x, h etc are numbers. 
But when we go to the surface of a balloon, we do not know how to implement such 
definitions eg how to take the 'difference of two points' on the surface. The way the 
idea of differentiation is captured is to assign numbers to points and then use them 
in the usual manner. This assignment can be thought of as pasting small pieces of a 
graph paper on the surface, and reading off the numbers corresponding to the points 
just below the numbers. This pasting gives a set of local coordinates in that portion 
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of the surface. Clearly there is huge freedom in the choice of the pieces of graph 
paper as well as in the pasting. Thus, while local coordinates can be introduced, the 
arbitrariness must be respected. Changing these assignments is called local coordinate 
transformations. Likewise, quantities such as tiny arrows stuck on the surface can be 
described in terms of components of the arrows with respect to the coordinate axes 
provided by the graph paper, When the graph paper is changed, these components 
change, but not the arrows themselves. Consequently, the two sets of components 
must be related to the local coordinate transformations in a specific manner so that 
both the sets refer to the same arrow. The arrows are an example of a vector field 
which can be thought of as a collection of their components, transforming in a specific 
manner. The components are denoted as quantities with arbitrary number of upper 
and lower indices. These must transform as, 

Objects represented by such indexed quantities have a meaning independent of the 
choice of local coordinates and are called Tensors of contravariant rank m and covari- 
ant rank n. Equations involving tensors (correctly matched index distribution) are 
covariant (or form invariant) with respect to local coordinate transformations. This is 
all that we will need to know. 

2. In the context of Lagrangian framework, the generalized coordinates will be denoted 
by q\ These are not necessarily Cartesian and are to be thought of as arbitrary local 
coordinates on the configuration space manifold. Consequently, we will occasionally 
also comment on whether various equations/conditions are 'covariant /invariant' under 
coordinate transformations. For these purposes, tensor notation will be used freely. La- 
grangian will always be taken to be a function of generalized coordinates and velocities 
and independent of time. 

3. In the context of Hamiltonian framework, the generalized canonical coordinates will 
be denoted as q\pi- It will be convenient to subsume them as 2n coordinates uA The 
phase space will be denoted as T which is a 2n dimensional manifold. On this there is 
a distinguished antisymmetric rank 2 tensor, VL^ v {uj) which satisfies further properties 
namely, 

dxVL^ + cyclic = 0. (1.3) 

The antisymmetric matrix is assumed to be invertible and its inverse is denoted 
as tt» x ttxu = 6g. 

Such an antisymmetric tensor is called a symplectic form. The coordinates in which 
the matrix takes the block off-diagonal form 



are the canonical coordinates and these are guaranteed to exist by Darboux theorem. 
4. The usual Poisson bracket of two functions f(cu),g(u) get expressed as, 

{f(uu),g(uj)} = ,f/ 9g 
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The ( II. 3p condition implies the Jacobi identity of Poisson brackets. 

For every function / on the phase space, one can associate the so-called Hamilto- 
nian vector field, Vj := VL^ u d u f . Such Hamiltonian vector fields generate infinitesimal 
canonical transformations as, 

u^-tv' 11 :=u» + evf(u). (1.6) 

It is easy to check that the symplectic form is invariant under these coordinate trans- 
formations (and hence these are called canonical transformations). In particular, the 
Hamiltonian function generates time dependent changes in coordinates which is noth- 
ing but the dynamical evolution. Hence, dynamical evolution can be viewed as a 
"continuous unfolding of canonical transformations". Taking e := 5t, f = H implies 
the Hamilton's equations of motion, 

= a-g. (i.7) 

dt duo" y ! 

Note that these are not just any system of first order, ordinary differential equation, 
but have a specific form involving the antisymmetric tensor VL^ U because of which H 
is constant along solutions of equations of motion. Furthermore, one also has the 
Liouville theorem regarding conservation of phase space volumes. 

We will begin our discussion by taking the example of source free electrodynamics. 
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Chapter 2 

Hamiltonian Formulation of 
Electrodynamics 



We will begin with the usual Maxwell equations, put them in the four dimensional rela- 
tivistic form, arrive at an action formulation from which we will go to the Lagrangian and 
Hamiltonian form. To be able to use the four dimensional relativistic tensor notation, we 
need to choose a set of conventions. 

Maxwell Equations: 

(2.1) 

3 (2-2) 

(2.3) 



V -E = p , Vx£ + 

V -B = , VxB- 
dp * -> 

— = -V • 7 

dt J 



dB 
~8t 
8E 
~dt 



Conventions: 



• Write all quantities without derivatives and with arrows as quantities with upper in- 
dices i.e. as contravariant tensor, eg, a B- (a 1 , a 2 , a 3 ) -B a 1 and vector derivatives as 

— * 

derivatives with lower indices (covariant tensor), eg, V B- (81,82,83) B- 8{. 

Thus, electric field, magnetic field and vector potential are contravariant tensors while 
the gradient operator is a covariant quantity. The cross product explicitly gives a 
contravariant quantity. 

• Quantities with upper indices are related to those with lower indices by, a 1 = — Oj, a = 
a i.e. a M := rf v a v with rf v = diag(l, —1, —1, —1) = rj^. 

• Cross products are expressed as, (V x a) 1 := e\dja k , where, e 123 = 1 = —6123. It then 
follows that, 

^^„ = -(M-^) • (2.4) 

Maxwell equations now become, 

BE 1 

8 % W = p , e\d,B k - — - = f (2.5) 
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diE* = , e ij k d j E k + — = (2.6) 

Define, j° := p, F 0i := —E\ F ij := -i\B k . It is easy to check that 

(1231) -> = , (2.7) 

(jZSD -»■ d>*F l/X + cyclic = = d^F u \ + cyclic (2.8) 

It follows that I?* = \^ % j^F^ k and the identification of the electric field matches with the 
usual definition with N 1 -H- (A := $, A). 

This allows us to think of Maxwell equations as tensor equations involving 4 dimensional 
tensors. If one transforms the coordinates x M as x' M := A^x u and transforms the F,j quanti- 
ties by tensor rules (index-by- index action), then evidently the Maxwell equations are form 
invariant or covariant for all invertible matrices A(J. Remember though that we have also 
used specific rules to relate the upper and the lower indexed quantities. These rules must 
also be respected by the primed observer (coordinates) i.e. rj^ must also be invariant under 
the coordinate transformations and this restricts the to be the familiar Lorentz transfor- 
mation matrices. We thus also see the Lorentz covariance of Maxwell equations. Incidentally, 
operating by <9 M on (I2.8P and using (12 .7p one can deduce the wave equation, DF^ U = and 
invariance of □ := r]^ v d^d v is precisely the requirement of Lorentz invariance. 

Note that the 4-tensor notation, is just a compact notation to write Maxwell equations and 
the physical property of electrodynamics being Lorentz covariant is encoded by the invariance 
of □ which in turn requires invariance of i]. The compact notation helps to keep track of the 
Lorentz covariance property. 

Now we would like to see if these equations can be obtained from an action principle and we 
all know that the answer is of course yes. For our purposes, it is sufficient to consider the 
source free case and we will now take = 0. 

If F^y are treated as basic variables describing electrodynamics, then we have 8 first order 
partial differential equations for six quantities. Usual equations of Lagrangian framework 
are second order (in t) equations for configuration space variables. However one notices that 
(12.81) can be identically solved by putting F^ u := d^A u — d u A^. This is always possible to do 
locally. Now the remaining 4 equations become 4 second order equations for the 4 quantities 
Ap. At least now the number of equations equals the number of unknowns. However, the 
definition of F in terms of A does not determine A uniquely (not even up to a constant); 
A'^ = A fl + <9 M A gives the same F. Therefore, although we have correct number of equations, 
these equations do not suffice to determine the candidate dynamical variables, A, uniquely. 
Modulo this observation, let us go ahead any way by thinking of ^(i") as 'configurations 
space variables'. 

It is a very easy exercise to check that if we define an action, 

S[A(x)} := JdtJ d 3 x (—F^fA := j d 4 x£(A„d a Ap) , (2.9) 

then its stationarity condition gives the Maxwell equations expressed in terms of A. This of 
course is an example of a field system i.e. a dynamical system with infinitely many degrees 
of freedom, A^t^x 1 ) with x l serving as labels. To proceed further, let us introduce some 
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working rules and notations: derivatives of the Lagrangian density etc will be denoted as 



SC(x) 
6A„(x>) 



with the rules, 

5A v {xf) 



5d u A /I (x) 
5A x (x>) 



SfoJ^x* - x*) 



(2.10) 



Now let us try to get to a Hamiltonian formulation. Our basic configuration space variables 
are A^t^x 1 ) while the Lagrangian is L = J d 3 x(— 1/4) F^F^. The conjugate momenta are 
given by 



:= / dV- 



5C{x) 



F 0fX (x) . 



(2.11) 



SdoA^x') 

Clearly tt° = identically while vr i = -F 0i . Therefore, m = -F 0i = -d Ai + diA . This 
allows us to express the 'velocities' d Ai = — 7Tj + diA which is needed in getting to the 
Hamiltonian form. However, we cannot do so for d Q A Q \ 



The canonical Hamiltonian is defined as, 

1 



d 3 x 
d 3 x 
d 3 x 



7r»d A„ + -F, V F^ 



(2.12) 



7T 



1 (doAi - diA + d t A ) - l - (F°f + X -F l3 F^ + ir°d A c 



(2.13) 



In the last equation, a partial integration has been done. The first two terms are the 
usual electromagnetic field energy density, E ■ E + B ■ B. The last term depends only on 
7T°, doA and is not a function of phase space variables, it contains a velocity. If we naively 
consider the Hamilton's equation of motion for A , we get an identity, implying that the 
time dependence of Aq is undetermined. The equation of motion for 7r° however leads to 
Oqito = diTT 1 . Observe that if we drop the last term by setting tt° = at some initial 
time, and require this condition to hold for all times, then we will need diir 1 = for all time 
which is just one of the Maxwell equations (the Gauss Law). This is also suggested by the 
Legendre transformation involved in going from the Lagrangian to the Hamiltonian - we are 
required to use the definition of the momenta. Hence, in the Hamiltonian we should not 
have the last term. Since tcq = from definition at all times, we must also obtain <9o7To = 
from the Hamilton's equation of motion. Thus it is self consistent to interpret the dropping 
of the last term as imposition of a constraint, tc° = on the phase space of (A^, 7r M ) and 
the Gauss law following as a consistency condition. Since A remains undetermined, we will 
replace A in the third term by an arbitrary function of time, X(t). Since the Gauss law is 
also independent of positions/velocities, we will interpret it also as a constraint, x := di 71 ' ''■ 
We will refer to tt° = as a primary constraint and x — as a secondary constraint, because 
it is needed for the primary constraint to hold for all times. 



Let us consider now the canonical transformations generated by the two constraints, 
follows, 



It 



SeA^x) 



d 3 x'e(x')n°{x') 



5 v A^x) = l^A^x), J d 3 x'r ] (x')d i Tr i (x' 
5n(x) = 



s°Ax) 

■ - S\dir)(x) 



(2.14) 

(2.15) 
(2.16) 
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Clearly, S e A (x) = e(x), 5 v Ai(x) = —di](x). Choosing e(x) := — dor)(x), reveals that the 
canonical transformations generated by the constraints are nothing but the usual gauge 
transformations. 

Let us summarize. 

1. Maxwell equations can be written in a manifestly relativistic form. 

2. These can be derived from relativistic action principle treating A^{t, x l ) as generalized 
coordinates. In the Lagrangian formulation, the 'matrix of second derivatives', 



is non-invertible. 

3. In the Hamiltonian formulation, there are constraints. These generate canonical trans- 
formations which are the familiar gauge transformations. 

4. The canonical Hamiltonian has the form which contains the constraints with arbitrary 
functions of time as coefficients. 

5. We also know that while manifest special relativistic formulation requires us to use 
4 components of the vector potential as basic configuration space variables (i.e. 4 
degrees of freedom per point), physically there are only 2 degrees of freedom (the 
two polarizations). Thus the constraints inferred above have something to do with 
identifying physical degrees of freedom. 

We will keep these in mind and try to abstract these features in a general framework. On 
the one hand we will simplify by restricting to finitely many degrees of freedom and at the 
same time generalize to arbitrarily complicated systems. 



5 2 C 



M^(x, x') = (-77^ + ?/V") S 3 (x - x') , (2.17) 
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Chapter 3 



Constrained Lagrangian and 
Hamiltonian Systems 



Consider a dynamical system with finitely many degrees of freedom, described in a La- 
grangian framework. Let Q be an n-dimensional manifold with local coordinates q l serving 
as generalized coordinates. Let L(q, q) be the Lagrangian function of the dynamical sys- 
tem and S[q(t)] := / dtL(q,q) being the action. The Variational Principle leads to the 
Euler-Lagrange equations of motion, 

d fdL\ dL . 1 n , n „. 

Tt\W) = 5? . = i,2,-,». (3.1) 

WW) * := ? = 5? - 5^ = «<<*•« ■ < 3 ' 2 ' 



Usually, the matrix M^- is invertible which allows us to solve for the accelerations, cf = 
(M- l ) ij Qj. This implies that given position and velocity at an instance, one can always 
determine the dynamical trajectory at other instances. 

If on the other hand the matrix is not invertible, then let 1 < r < n, denote its 
rank. Then there exist (n — r) independent h l a , vectors satisfying Mh = and we cannot 
solve for all the accelerations uniquely. The general solution for the accelerations will be, 
q l = q l + Yla=i a oh % a an d the first term is the solution of the inhomogeneous equation. 
Furthermore, Mh = implies h T M = = h l a Qi(q, q) — for a — 1, • • • (n — r). Thus on the 
one hand the equations of motion are not uniquely specified and on the other hand there are 
(n — r) relations among the 2n positions and velocities. If these relations are not satisfied 
identically, then in particular, they correspond to restrictions on the initial data, g l (0), q l (0). 

As an example, consider a system with a one dimensional configuration space with coordinate 
q and with a Lagrangian, L = f{q)q — V(q). Clearly, M = and Q = —V'(q). The 
acceleration is given by q = ah and hQ = =>• V'(q) = 0. If V = then the equation is 
satisfied identically while for non-zero V, the position must be at an extremum of V (which 
may not exist!). 

Let us see what implications are there for a corresponding Hamiltonian framework. First 
observation is that for a singular Lagrangian, the definition of conjugate momentum p { := ^ 
cannot be inverted to solve for the velocities in terms of the momenta. This is essentially the 
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inverse function theorem namely, y % = p(x) can be inverted to get x % = g z (y) provided is 
an invertible matrix. (Strictly speaking, the n momenta are n functions of the 2n positions 
and velocities and hence one should be invoking the implicit function theorem, such fine 
prints are ignored here.) 

Let us define the Hamiltonian by H := piq % — L(q, q) and consider its variation, 

Pi ~ W ) & - Q^M ( 3 - 3 ) 

If we treat the q, q, p all as independent variables, then the 'Hamiltonian' is clearly a function 
of all of these. Let us use the definition of momenta (i.e. restrict to a 2n dimensional sub- 
manifold of the 3n dimensional space). On this sub-manifold, the middle term vanishes and 
the Hamiltonian becomes a function only of positions and momenta in the sense that it varies 
only when q i ,pi are varied. This is independent of whether velocities can be solved for in 
terms of momenta and positions. The next step in the usual case of non-singular Lagrangians 
is when one infers the Hamilton's equations of motion as, 

q* ~ dH v ~ - dH (3 4 ) 

dpi ' 1 dq l 

This assumes that (a) the variations 5pi, 5q l are independent and (b) Hamiltonian is a function 
only of Pi,q l . Both these statements fail for singular Lagrangians. 

To see this, notice that variation of the Hamiltonian being given in terms of variations of 
momenta and positions, depends on evaluating the variation of the Hamiltonian on the sub- 
manifold defined by Pi = |tj. If the variations are also to respect this condition, we must 
have 

Sp, = M v Sgi + -^-S q j (3.5) 

which immediately shows that the variations Spi, 8q l are not independent. Consequently, one 
cannot deduce the Hamilton's equations of motion. We have already noted that for singular 
Lagrangians, the velocities cannot be eliminated in favour of momenta and positions. The 
variations being not independent means that there are relations among the phase space 
variables Pi,q\ namely, (p a {q,p) = 0. 

Let us summarize. 

A Lagrangian system is said to be singular if the matrix of second derivatives of the 
Lagrangian with respect to the velocities, is singular. This has the consequences that (a) the 
accelerations are not determined uniquely i.e. contain arbitrary functions of time and (b) 
it could imply relations among velocities and positions which may not even be consistent. 
It is possible to isolate the accelerations which are determined uniquely. However we will 
carry out such analysis in the context of a Hamiltonian framework. The property of being a 
singular Lagrangian is independent of any choice of local coordinates. Generically, singular 
Lagrangian systems are also called constrained systems 

In a Hamiltonian formulation obtained from a Lagrangian, the singular nature of the La- 
grangian manifests as certain relations among the phase space variables. We will therefore 
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define a constrained Hamiltonian System to be a system with a 2n dimensional phase space 
manifold T with local coordinates {q l ,Pi) -H- on which is given a Hamiltonian function, 
H (q l ,pi) := H (u) together with a set of relations <p a {u) = 0, a = 1, • • • k(< 2n). These rela- 
tions are referred to as Primary Constraints. By definition, the k constraints are functionally 
independent i.e. the k differentials, d(j) a (id) (or the k vectors, d^(j) a ) are linearly independent. 
The sub-manifold defined by <p a = will be denoted by S and is called the constraint surface. 
It is a In — dimensional sub-manifold of T. Note that a constraint surface is not a phase 
space in general, i.e. does not have a symplectic form (eg when k is odd). 

We will focus entirely on the constrained Hamiltonian systems and analyze various possibil- 
ities of types of such systems. The aim will be to have a procedure for obtaining equations 
of motion in the Hamiltonian form paying attention to the constraints. This means that we 
want to have a variational principle for paths in T, which will lead to 'Hamilton's equations 
of motion' for some suitable Hamiltonian function and such that the possible dynamical tra- 
jectories (in T) either remain confined to the constraint surface, S or never intersect it. In 
effect, we can continue to work in the given phase space and use a new Hamiltonian function 
so that dynamics effectively respects the constraints. This is achieved by using the method 
of Lagrange multipliers. 



S[u(t),X(t)] := / dt 



-y% u u v - H (cu) - \ a {t)<p a {u) 



(3.7) 



Introduce k Lagrange multipliers, A a and define a new Hamiltonian function H := H + \ a (p a , 
which matches with the given Hamiltonian H on the constraint surface, S. Defining an 
action, 

"1 

2 

and invoking its stationarity, 5S = 0, leads to the equations of motion, 

^^(^ + ^) md *>) = <>, .= 1,-,*. (3-8) 

Thus we obtain equations of motion in a Hamiltonian form and also the constraint equations. 
We have succeeded in having a new Hamiltonian dynamics defined for trajectories in T. We 
have now to ensure that the trajectories are such that if an initial point is on the constrained 
surface, then the whole trajectory also remain on the constrained surface. If this property 
can be ensured, then it also follows that no trajectory can enter and leave E, since the 
equations are first order. 

Observe that if the value of <p a for any given a is preserved under evolution (i.e. <p a = 
along a trajectory), then the trajectory is confined to the 2n — 1 dimensional hyper-surface 
cf) a =constant. However we only need the trajectory to be confined to S, so the value of 
each constraint need not be preserved exactly. Therefore we need not have <p a = along all 
trajectories but only along those trajectories which lie in S. This is ensured by requiring 
that the Poisson bracket of the constraints with the Hamiltonian H be weakly zero. Since 
Poisson brackets can be evaluated at any point of the phase space, we can evaluate these at 
points on the constraint surface and weak equality/equations refer to Poisson brackets being 
evaluated at points of E and are denoted by '~' [2]. Strong equations/equalities are valid in 
a neighbourhood of £ and are denoted by the usual '='. In particular a strongly vanishing 
function vanishes weakly and so do all its partial derivatives. 

^ = {r,H } + \ b {<p a ,<P b } « (3.9) 

These k conditions ensure that a trajectory beginning on the constrained surface remains on 
the constraint surface and our goal is reached provided (13.91) holds on S. 
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There are several possibilities now [2]. 

If the k x k matrix of Poisson brackets of the constraints is non-vanishing, then the con- 
sistency conditions can be viewed as a matrix equation at each point on X for the Lagrange 
multipliers X a . Since this matrix is antisymmetric, it is non-singular only if when k is even. 
In this case, all Lagrange multipliers are determined and we do have a Hamiltonian dynam- 
ics whose trajectories either lie on the constraint surface or never intersect it. Generically 
however the matrix is singular. As seen in the context of singular Lagrangian, this means 
that (a) some multipliers are necessarily undetermined and (b) some linear combinations of 
{0 a ,ifo} must vanish on E. Again we have several possibilities. Either (i) {cf) a ,Ho} « 0,Va 
and all linear combinations are weakly zero, or (ii) the linear combinations vanish weakly 
provided some further functions vanish in which case we refer to these as secondary con- 
straints, or (iii) there are no points of S at which the linear combinations vanish in which 
case we say that the Hamiltonian system is inconsistent. 

It could also happen that the matrix is zero on the constraint surface. This could happen 
for instance, if {4> a , <j) b } = C ab c cf) c . (Recall the Maxwell example). In such a case, there is no 
equation for the Lagrange multipliers and all Lagrange multipliers are undetermined. This 
case can be thought of as a special case of rank of the matrix being zero. 

If we encounter the inconsistent case, we throw away our formulation of the system and 
start all over again. In the case (i), we have reached our goal but have to live with some 
undetermined Lagrange multipliers (and hence evolution). This will turn out to be the most 
interesting case. In the case (ii), we have to now demand that the Poisson bracket of sec- 
ondary constraints with the Hamiltonian H must vanish on S. Once again we will encounter 
similar above and we have to repeat the analysis - we will either satisfy the condi- 

tions identically with some Lagrange multipliers determined or encounter tertiary constraints 
or encounter inconsistency. Since the total number of constraints cannot be larger that 2n 
(else no initial condition will be left!), the process must terminate. Barring inconsistent 
systems, we will eventually end up with some Lagrange multipliers being determined, some 
undetermined and with possibly additional constraints xa ~ 0, A — 1, ■ ■ ■ , I. 

Note: We began by requiring the trajectories to be confined to £ and found as a consistency 
requirement that the goal cannot be attained for all points o/S. We need to restrict further 
to a sub-manifold £' C S, due to the secondary constraints. Thus, the dynamics defined by 
H on V , will correspond to a constrained dynamics relative to £' defined by all constraints 
being weakly zero. The dimension of £' is of course 2n — k — l. Notice that even if we began by 
requiring consistency condition to hold on E, extending it to hold on does not contradict 
the previous condition since £' is a sub- manifold of S. The weak/strong equation now refer 
to E'. 

To summarize: Beginning with a phase space T, a Hamiltonian Hq and a set of primary 
constraints a , we can construct a Hamiltonian dynamics on T such that its trajectories 
are either confined to the constraint surface £ or avoid it. The construction reveals the 
possibility of further constraints as well as the dynamics being not completely determined. 
In the next lecture, we will consider a suitable classification of constraints and obtain a 
corresponding classification of constrained Hamiltonian systems. 
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Here are some elementary examples. 

2 

1. H = 7^, (f)(q) = q — g : {(f), H} = ^ ^ , p ~ 0, is a secondary constraint; 

2. H Q = fi 0(g) = g 2 : {0, if} = ^ » condition holds; 

3. ffo = ap, a 7^ 0, 0(g) = g : {0, ff } = a ps inconsistency. 
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Chapter 4 



DiraoBergmann theory of 
Constrained Hamiltonian Systems 

Our consistent constrained Hamiltonian system is specified by, 

k 

H = # o + ^A a a , a ^O , x A ~0, A=l,---l, k + l<2n. (4.1) 

a=l 

where £ ~' means evaluation on X defined by the primary constraints (<j) a ~ 0) and the 
secondary constraints (x A ~ 0). The constraints also have to satisfy, 

{0 a ,/J o } + {0 a ,0 fe }A b « , {x A ,H } + { X A A b }h ~ 0. (4.2) 

Thus we have + / equations for Lagrange multipliers and the system is naively, over- 
determined. 

Let rank of the (k + I) x k matrix of Poisson brackets of the constraints be K < k. This 
means that K is the maximum number of linearly independent rows and columns of the 
matrix. Thus there exist k — K independent relations among the k columns of the matrix 
i.e. 3 ^ b a \ a = 1, • • • , k — K numbers such that 

{0 a ,0 b Ki Q) ^O , { X A A b }& ] ^0. (4.3) 

Now define two sets of linear combinations of the primary constraints namely, 

r := , ¥ := Vi a ' ] r , a = l,---,k-K and a' = l r --,K. (4.4) 

The new set of K vectors i]i a ^ are so chosen that the the set of constraints <fi a , 4> a ' are 
functionally independent. Now it is clear that 

{r,r}^o , { x A A a }^o. (4.5) 

Thus the k — K new combinations, <f) a of primary constraints have a weakly vanishing 
Poisson bracket with all the constraints. Such constraints are termed first class constraints. 
Constraints which do not have this property are termed second class constraints. With the 
help of the £'s we have regrouped the primary constraints into primary first class and primary 
second class constraints. 
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We would like to do the same for secondary constraints. Observe that linear combination of 
second class constraints will again satisfy the consistency conditions and to such combinations 
could be added any combination of primary constraints without affecting these equations. 
To maintain the functional independence of the secondary constraints we need to ensure that 
the linear combinations are also functionally independent. Thus, consider the combinations, 

X A := S A BX B + SU a + SU a ' (4-6) 

with S' . so chosen that S B is a non-singular matrix and x A have weakly vanishing Poisson 
bracket with all constraints for a maximal number of values of A. Let this number be L < I. 
Thus, we divide the combinations x into the first class combinations, a = 1, • • • , L and 
the second class combinations c' = I, ■ ■ ■ ,1 — L. 

The result of these manipulations is that (a) we can write \ a (f) a = \ a 4> a + A<y0 a ' and (b) the 
consistency conditions can be simplified as, 

{r,H }^0 , {r,H }^0 (4.7) 

{r',H } + v{0 a ',^'}^o , {r',#o}+v{r',^'}~o. ( 4. 8) 

The first set of equations involve only the k — K + L first class constraints and no Lagrange 
multipliers while the second set of K + I — L equations involve only the K Lagrange mul- 
tipliers, \ a i. The k — K Lagrange multipliers, X a have dropped out of the equations and 
will remain undetermined. Once again we have more equations than unknown, but because 
of the separation into first and second class constraints, we are now guaranteed that the 
(K + 1 — L) x K matrix of Poisson brackets of (p a and with <^r has the maximal rank K . 
For, if it did not, there would exist further linear combinations which will weakly Poisson 
commute with all constraints and by construction we have obtained the maximum number 
of first class constraints. We will now solve for the Lagrange multipliers X a > explicitly. 



Define the matrix A of Poisson brackets of the second class constraints as, 

(4.9) 



{4> a ',^'} {4> a ',x a '}\ /A '/ 3 ' a qV 

{x p \r'}) : " 

This is an antisymmetric square matrix of order K + / — L. This must be non-singular. 
For, if it is singular, there will exist non-trivial linear combination of the second constraints 
(p a , x a , which will weakly Poisson commute with all the second class constraints (and it 
automatically commutes with the first class constraints), implying that we have additional 
first class constraint. The non- singularity also requires that the total number of second class 
constraints, K + I — L, must be an even integer. Let its (weak) inverse be the matrix C, 

C := ( C n a ' p ' ) (4.10) 



Cpi/3' C p i a i 



The equation CA w 1 translates into, 



C^A^ + CWA^ = V (4.11) 

C^A^' + C^A ctV = 5/ (4.12) 

C a , p ,^' T ' + CWA ctV = (4.13) 

C^A^ + CpvA^ = (4.14) 
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The equations for the Lagrange multipliers become, 

{^,ff o } + A^V»0 , {r',#o} + A^'A 7 ,«0 (4.15) 

Multiply the first one by C a ip>, second one by C a ' a >, add the two and use (14. lip to solve for 
\ a >. One gets, 

V « -C a ^{0 /3 ',i/o}-C a v{x' J ',^o} (4.16) 
Similar manipulation using (I4.14p equation leads to, 

C p , ,{4>P',H o } + C p/a ,{x a ',H o } « (4.17) 

We can now write the Hamiltonian explicitly using the solution (I4.16P and express the 
evolution of any function on the phase space as its Poisson bracket with H. We will write it 
in a more convenient form. 

j t fHt)) « {f,H } + X a {fA a } 

-{f,4> a '}C a ,p,{^\H } - {/,f'}(7 aV {f> } (4.18) 
-{f,X p '}C P >p>{^\H } - {f,r'}C P > a/ {r',H } 

The last line, which is weakly zero due to ( 14. 17ft . has been added to get a more symmetric 
final expression. 

Now let us denote all the second class constraints, (<fr a ' , x p ') by £ m , m = 1, • • • , K+l—L. Then, 
the nonsingular matrix A is just the matrix A mn = {£ m , £ ra }, C mn is its weak inverse as before 
and the last group of four terms in ( 14. 18ft are conveniently expressed as — {/, £, m }C mn {!; n , H } 
so that finally we obtain, 

« {/, ff } + ^ A a {/, f } - {/, rKA-^nir, i?o} , (4.19) 

a 

Q 

= {CO and (4.20) 

= {f,9}-{f,C}(^- 1 )mn{C,9} (Dirac Bracket) . (4.21) 

In this final expression, we have removed the ~, the primary first class constraints are denoted 
by <p a while all second class constraints are denoted by £ m . 

Several remarks are in order. 

1. The first step in the analysis of constrained systems was to obtain the full set of 
constraints starting with a given Hamiltonian Ho and a set of primary constraints 
defining the constraint surface S C T. In order to ensure that we get the final form of 
evolution equations to be a Hamiltonian form, we used a modified Hamiltonian H and 
thought of the system as thought it were un- constrained in the sense the variations of 
the phase space coordinates were independent. To make contact with the constrained 
nature of the system, we required that the un-constrained dynamics be such that its 
trajectories either lie in £ or never intersect it. This lead us to discovering possible, 
additional constraints. Note that restrictions on the trajectories was with reference to 



A mn 
{/,<?}* 
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the sub-manifold £ and hence only the primary constraints played a role in subsequent 
analysis i.e. we did not add to the Hamiltonian, terms corresponding to the secondary 
constraints. The secondary constraints however do reveal that the required segregation 
of trajectories holds only with respect £' C X, defined by vanishing of all constraints. 

2. The next step was essentially an exercise in linear algebra. We did this to solve ex- 
plicitly for those Lagrange multipliers which could be solved for. This was facilitated 
by regrouping the set of all constraints into first class and second class constraints. 
The final result reveals that evolution could be arbitrary if there are primary, first 
class constraints due to the undetermined X a . The most compact expression for the 
Hamiltonian evolution was obtained using the Dirac brackets. 

3. We now define first class variables as those functions on T whose Poisson brackets with 
all constraints are weakly zero. By the consistency condition (14. 2p . the Hamiltonian H 
is a first class variable and of course so are the first class constraints. The Hamiltonian 
H may or may not be a first class variable. It's Poisson bracket with first class 
constraints is of course weakly zero from (14. 7p . It is easy to check that sums and 
products of first class variables is again first class and so are the Poisson brackets of 
first class variables 0. 

4. As noted already, in the presence of primary first class constraints, evolution equation 
for a generic function /, contains the arbitrary Lagrange multipliers, \ a . From (I4.19p . 
it follows that evolution of first class variables is independent of X a and is entirely 
governed by H . This justifies why first class variables are singled out. 

5. The Dirac brackets have been introduced as a convenient compact notation. However 
it has many interesting properties, namely, 



Thus it has all the properties of the usual Poisson bracket. In addition, one has 

(a) {f,g}* ~ {f,g} for any first class / and \/g; 

(b) {C,g}* ~ 0, Wg and any second class constraint £\ 

Recall that weak equations involving Poisson brackets mean that the Poisson brackets 
are first computed in a neighbourhood of £ and then evaluated on E. This rule is 
necessary since a weakly vanishing function need not have a weakly vanishing Poisson 
bracket (since some of the partial derivatives 'off' £ may not be zero). This applies to 
the second class constraints as well. However, Dirac bracket of a second class constraint 
with any function is weakly zero. Therefore, if we use Dirac brackets for writing the 
equations of motion (as shown in (I4.19p ). then we can set the second class constraints to 
be zero before computing the Dirac brackets. This is equivalent to reducing the phase 

functions on the phase space of an unconstrained system are generally called observables while in the 
context of a gauge system, the first class observables are all called as Dirac observables. 




{f,9 + h}* 
{f,M}* « 

{Lgh}* « 



{f,g}* « 

{f,{9,m 



{f,g}* + {f,h}* (addition); 
[i{f,g}* (scalar multiplication); 
{f,g}*h + g{f,h}* (Leibniz); 
- {gJ}* (antisymmetry); 
+ cyclic ps (Jacobi identity). 
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space dimension from 2n to (2n - the number of second class constraints). Thus, second 
class constraints correspond to redundant degrees of freedom which can be ignored by 
setting them to zero. 

6. We could now focus on systems that do not have any second class constraints either a 
priori or after eliminating them via the Dirac bracket procedure. One is effectively left 
systems with only first class constraints. As noted earlier, the evolution of a dynamical 
variable in such systems is in general, arbitrary and only variables whose evolution is 
not arbitrary are the first class variables for whom the Dirac brackets are same as the 
Poisson brackets. 

Note that it could happen that there are no first class constraints left. In such a case, 
we have just the usual types of systems. Thus, the net conclusion of the analysis is: 

Generically, consistent Hamiltonian systems are systems with first class constraints 
with the special case of no first class constraints. Hamiltonian systems with at least one 
first class constraint are termed gauge theories. We will now focus on these exclusively. 

Let r be a 2n dimensional phase on which is given a first class Hamiltonian function Hq and 
a set of of k < 2n first class constraints, <p a whose vanishing defines the constraint surface 
E. The total Hamiltonian governing time evolutions is given by H := H + J2 a ^ a (f )a - By 
virtue of being first class, we have the following relations: 

{0 a , 6 } « <-> {<P\ <f) b } = C ab c (uj)<p c and {H , a } « <-> {H , <p a } = D\(uj)<p b . 

(4.22) 

The evolution equation for any function / on Y is given by, 

j t f{u{t)) « {f(co),H + J2^r(co)}l = ^t) (4.23) 

a 

Now observe that (a) if we make an arbitrary diffeomorphism i.e. a mapping of the manifold 
r on to itself preserving differential structure, the manifold is unchanged (by definition). 
However, the symplectic form would change in general; (b) if we specialize the diffeomor- 
phisms to those which preserve the symplectic form, then the restricted diffeomorphisms are 
the familiar canonical transformations. All such (continuously connected to identity) trans- 
formations can be generated by arbitrary functions on T, by the rule: 5 eg u^ := eWd v g{u). 
All of these however do not preserve the constraint surface; (c) To preserve the constraint 
surface, the function must preserve the constraints defining the surface i.e. must be a first 
class function. Thus, all first class functions and in particular the first class constraints, do 
preserve E. While constraints preserve the Hamiltonian, other first class functions need not. 
Those first class functions which do preserve the Hamiltonian as well are said to generate 
symmetry transformations. By contrast, the transformations generated by the first class 
constraints are distinguished as gauge transformations. 

Thus, the diffeomorphisms generated by first class constraints, preserve the entire structure of 
the constrained Hamiltonian system (i.e. manifold, symplectic structure, constraint surface 
and the Hamiltonian) and are termed gauge transformations. First class functions (which 
are not constraints) are automatically invariant under these transformations. Provided they 
Poisson commute with the Hamiltonian, they generate symmetry transformations. This 
distinction among the set of all first class function, comes about for the following reason. 
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Recall that A a are undetermined, arbitrary functions that appear in the equations of motions. 
Therefore, beginning from any initial condition Co G £, the actual trajectories will depend 
on the choice made for the A a 's. //we identified points on E as representing physical states 
of the system, we would loose determinism - a given state does not uniquely determine the 
future state. We need to identify physical states of the system differently. 

Consider infinitesimal evolutions for two different choices of A a 's. We will have, 

u'(6t) = u + 6t{u,H(X',uj)}\c , u(5t) = u) + 5t{u,H(\,u)}\a , (4.24) 
which implies that, 

6u(6t) = 5t{u,J2^ar}\* = 5>^A a ){u;,0 a }U , (4.25) 

a a 

which is nothing but the infinitesimal transformation generated by the first class constraints! 

Thus, if we arbitrarily choose the A a and consider a trajectory evolved from some Co then 
another trajectory evolved by a different choice of A a from the same initial point, would be 
obtained by making a gauge transformation. Clearly, if we identified points in S which are 
related by gauge transformations as being 'physically the same', then we regain determin- 
ism in the sense that physical states evolve into unique physical states. Thus the apparent 
dynamical indeterminism implied by the first class constraints appearing in the Hamiltonian, 
can be resolved by defining equivalence classes of points o/E under transformations generated 
by the constraints (i.e. gauge transformations) to represent the physical states of the system. 
Notice that this identification of physical states with equivalence classes under gauge trans- 
formations involves only those first class constraints which appear in the Hamiltonian since 
only these have a bearing on the dynamical evolution. 

Since there are k first class constraints, the set of points of £ which are related by gauge 
transformations is parameterized by k parameters and hence the set of gauge equivalence 
classes is parameterized by 2n — k — k parameters. This space is called the Reduced Phase 
Space. This space be made explicit by introducing additional k 'constraints' (now called 
as gauge fixing conditions - Xa(w) ~ 0). These are required to be such that the matrix of 
Poisson brackets of the <p a , Xb constraints is non-singular. Demanding their preservation in 
time fixes the Lagrange multipliers and hence the name. 

Having clarified the identification of physical states, the definition of physical observables 
and notions of symmetry obviously must be formulated for physical states. The observables 
must have unique evolutions since by definition these are supposed to be functions of physical 
states. Only first class functions satisfy this property and only these qualify to be termed as 
physical observables. Likewise, notion of symmetry should refer to invariance with respect 
to transformations of physical states, the generators of infinitesimal symmetries must map 
the entire gauge equivalence classes among themselves and of course preserve the evolution. 
Clearly these again have to be first class functions and must additionally Poisson commute 
with the Hamiltonian. 

To summarize: 

1. Hamiltonian systems with at least one first class constraint, require identification of 
physical states not with individual points of the constraint surface but with gauge equiv- 
alence classes of points of the constraint surface. 
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2. In view of the above, it is common to refer to the original phase space as the kinematical 
phase space, T kin . The constraint sub-manifold of r kin is E. The physical state space 
(or reduced phase space) is denoted as r phys := E/ ~ where ~ refers to the gauge 
transformations. Note that the physical state space is not a sub-manifold of E. 

3. Although both the constraints and 'conserved quantities' serve to confine the tra- 
jectories the two are distinguished by the fact that constraint impose limitation on 
possible initial conditions (restriction to E) as well as force a non-trivial identification 
of physical states to ensure determinism of dynamics. Notions of conserved quantities, 
symmetries become meaningful only after this identification. Also conserved quantities 
do not impose any ah initio limitation on the possible initial conditions but only on a 
subsequent trajectory. 

Observe that in the light of the discussion of symmetry and gauge transformations, a first 
class Hamiltonian Hq generates a symmetry transformation, namely time translations. How- 
ever there are theories in which Ho = and H is made up entirely of first class constraints. 
Now the 'time evolution' itself becomes a gauge transformation and hence 'no evolution' of 
physical states. Next chapter discusses this case. 
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Chapter 5 

Systems with the Hamiltonian as a 
constraint 



Consider special types of gauge theories in which Hq = i.e. the Hamiltonian is entirely made 
up of first class constraints. The prime physical example of such a system is the dynamics of 
inhomogeneous cosmological space-times within the context of Einstein's theory of General 
Relativity 0. 

As is well known, Einstein's general relativistic theory of gravity (GR) has a four dimensional 
manifold on which is defined some metric tensor (of Minkowskian signature) field which 
makes it in to a space-time. The metric tensor is not a fixed entity, as in the case of 
special relativity (the Minkowski space-time), but is a dynamical entity i.e. is determined 
in conjunction with the (interacting) matter distribution on the manifold. The equation 
determining the metric is the Einstein equation, which is a set of 10, local, partial differential 
equations of order 2, for the 10 components of the metric tensor, g^ u . It is non-trivial fact 
that these equations admit a well-posed initial value formulation i.e. (i) one can take the 4 
dimensional manifold as IR x S3, (ii) specify two symmetric tensor fields, gij, and its time 
derivative, on S3, then the space-time can be determined for other times provided the 
'initial data' satisfies certain conditions. Furthermore, the system of equations can be cast in 
the form of constrained Hamiltonian system, with Hamiltonian given entirely by first class 
constraints. There are 4 sets of constraints (per point, since one is dealing with a field theory), 
three of which, called vector or diffeomorphism constraints and the remaining one called the 
scalar or Hamiltonian constraint. The vector constraints generate usual diffeomorphisms 
of S 3 while the scalar constraint generates evolution of the 'spatial manifold', S 3 in the 
4-manifold constructing a solution (space-time) of the Einstein equation. The upshot is 
that solution space-times of Einstein equation can be viewed as phase space Hamiltonian 
trajectories in the gravitational phase space with initial data satisfying a set of first class 
constraints and with the Hamiltonian given as linear combination of these constraints. The 
constraints in the Hamiltonian formulation of general relativity, reflect the 4-diffeomorphism 
invariance of GR [I]. 

This is pretty complicated to deal with in general, however a simplification is possible. If 
we restrict ourselves do the dynamics of only a special class of 3-geometries, namely, those 
metric tensors whose dependence on spatial coordinates (coordinates on S 3 ) is completely 



1 For space-times corresponding to compact objects, typically asymptotically flat space times, there is a 
'true Hamiltonian', the analogue of Hq, generating asymptotic time translations. 
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fixed in a particular manner and only time dependence is to determined, then GR simplifies 
to a Hamiltonian system for a finitely many degrees of freedom (6 for general homogeneous 
geometries, 3 for so-called diagonalized models and 1 for homogeneous, isotropic geometry 
of the Friedmann- Robert son- Walker (FRW) cosmology) with only the single Hamiltonian 
constraint remaining. We will not need to take any specific model to illustrate the issues. 

We can also construct systems in which Hamiltonian is the single (and hence first class) 
constraint. To see this, let us begin with a usual un- constrained Hamiltonian system with a 
phase space T and a Hamiltonian H (u). Let us extend this phase space to T by adding two 
more conjugate variables, r, tt. Let </>(r, tt, uS) := Hq(u) — 7r be chosen as the Hamiltonian on 
T and take it as a constraint as well, i.e. H := \cf) w 0. Clearly, any evolution (with respect 
to T) generated by the Hamiltonian is a gauge transformation and therefore unphysical. 
Functions which are insensitive to the evolution are the first class ones (Dirac observables) , 
/. Let us look for Dirac observables. 

{f,H} « {f(r,7T,u),H (uj)}-^ « 0. (5.1) 

where we have used the definition of Poisson bracket in the second equation. 

It is clear that functions which are independent of r, tt and satisfying the usual un-constrained 
evolution equation in T, are Dirac observables of the extended constrained dynamics. Func- 
tions depending only on r are not Dirac observables while those depending only on tt are 
Dirac observables. Constants with respect to the r-dynamics, also are Dirac observables. 

With this construction, we see that usual unconstrained dynamics can be viewed (albeit 
trivially) as a constrained dynamics in a bigger phase space. Furthermore, all functions on the 
unconstrained phase space, evolving by the un-constrained dynamics are Dirac observables 
of the constrained dynamics. The Dirac observables however are constants with respect to 
the T-evolution. 

Consider now the trajectories of the constrained dynamics. One finds, 

aW T = ~ X ' IT* = ° ' 5f W = M",#o(")} = A^c. (5.2) 

The last equality is deduced from ( 15. II) with f = u. We can define a new 'time', T" by the 
equation dT' = XdT, so that 

d d d d 

dT T = - 1 ' ^ = ° ' dT' U = fr U ' (5 ' 3) 

Note that T evolution is generated by \<p while T' evolution is generated by <fi. Starting 
with some initial values at T' = 0, one generates the T'-trajectories. All the points along 
these are related by gauge transformations generated by 0. Thus the equivalence classes 
are in one-to-one and on-to correspondence with the points r = 0, tt = tt, u> = Cj. The 
gauge orbits which lie on the constraint surface satisfy tt = H(u) and hence these orbits are 
completely determined by points in T. Thus the reduced phase space in this case is just the 
un-constrained phase space T. 

This simple construction brings out a few points. There are two notions of 'evolution' (i) 
the T or (T") evolution, which is some times called an external time evolution and (ii) the r 
evolution which is correspondingly called an internal time evolution. From the point of view 
of the constrained system, r is just one of the degrees of freedom which is singled out because 
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the constraint had a particularly simple additive form. The r evolution can thus be thought 
of as evolution of a set of degrees of freedom with respect to a singled out degree of freedom. 
The internal time evolution is thus also called a relational evolution while the singled out 
degree of freedom is called a clock degree of freedom. The arbitrary function, A is also called 
a lapse function and its arbitrariness corresponds to the freedom of re-parameterizing the 
external time. The Hamiltonian systems resulting from homogeneous cosmologies of GR, 
typically get presented in the form of the constrained system (the constraint however has a 
different form in general). The terminology used above is inherited from the GR context. One 
can in fact do a more general and systematic analysis of the notions of external and internal 
dynamics which my student Golam Hossain and I have carried out for finite dimensional 
systems. 

While classically constrained Hamiltonian systems are interesting in their own right, they 
become more challenging at the quantum level. As all the fundamental interactions of stan- 
dard model and its extensions are gauge (field) theories, one has to face these challenges. In 
the perturbative analysis, one needs to 'fix a gauge', in order that propagators for gauge fields 
can be defined and then has to show that the final observable scattering cross-sections are 
indeed gauge invariant. When a quantum theory of gravity is attempted, the understanding 
of semiclassical approximation becomes quite complicated especially in a non-perturbative 
approach. 
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Chapter 6 



Exercises 



1. Check that the coordinate transformations generated by Hamiltonian vector fields leave 
the symplectic form invariant. 

2. Check that ( II. 3p property is needed to prove the Jacobi identity for Poisson brackets. 

3. For the Maxwell theory, check that the secondary constraint holds for all times without 
having to require any further constraint. Furthermore, the Poisson bracket of the two 
constraints is also zero. 

4. Show that the Hamilton's equations of motion can be identified with the Maxwell 
equations. 

5. For the Maxwell theory, we have already shown that the canonical transformations gen- 
erated by the first class constraints are indeed the usual gauge transformations (hence 
in fact the name). Show, by direct computation, that are first class functions. 

6. Consider a massive, relativistic particle with action S = mo / dr ^Jrj IJLV x^x v . Carry out 
the constraint analysis and give examples of first class functions. 

7. Consider the four dimensional phase space with coordinates {q 1 ,q 2 ,Pi,P2)- Consider 
two constraints (f)(q,p) '■— p\ + P 2 , + (q 1 ) 2 + (<? 2 ) 2 — R 2 andx(?,p) '■— P2- Let H (q,p) be 
some suitable Hamiltonian (not given explicitly) such that these two constraints are 
preserved [5]. 

(a) Identify the constraint surface. 

(b) Compute the expression for the Dirac bracket. 
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